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摘 要 : 以 河 湖 为 主体 的 塔里木 河 下 游 地 表 水 体 的 盏 
水 平衡 中 扮演 着 重要 角色 。 基 于 近 20 a 塔里木 河 下 游 Landsat 5、Landsat 7 和 Landsat 8 影像 资料 和 
生态 输 水 数据 ,综合 采用 Google Earth Engine(GEE) 计 算 平台 和 多 元 统计 分 析 方 法 ,对 地 表 水 体 鱼 
变化 及 其 对 生态 输 水 的 响应 进行 了 系统 分 析 。 结 果 表 明 :(1) 2000 一 2019 年 ,塔里木 河 下 游 区 域 
积 均 呈现 波动 上 升 趋势 ,速率 分 别 为 15.54 km a, 


的 地 表 水 体 、 季 节 性 水 体 和 永久 性 水 体 的 


积 受 到 输 水 工程 的 影响 ,在 区 域 水 循环 、 调 节 


7.12 km?+a' 和 8.41 km”…a'; 塔 里 木 河 下 游 上 上 段 (大 西海 子 水 库 一 英 苏 ,不 包括 大 西海 子 水 库 )、 中 段 
( 英 苏 一 阿拉 干 ) 和 下 段 (阿拉 干 一 台 特 玛 湖 ) 区 域 的 季节 性 水 体 和 永久 水 体 


只 也 均 呈 现 波动 增 


加 趋势 ,其 中 下 段 区 域 的 增加 速率 最 大 ,其 值 分 别 为 5.23 km-a $8 824 km-a (2) 生态 输 水 引起 


BY He He ARE 


38 Jm Y 417.08 km'( # 10.92 È) ,(3) 近 20a 来 ,研究 区 域 地 表 水 体面 积 ,尤其 是 永久 水 体 


积 增加 主要 表现 在 台 特 玛 湖 湖 区 ,2019 年 台 特 玛 湖 湖 区 的 地 表 永 久 水 体 和 季节 性 水 
体面 积分 别 约 为 267.27 km 和 188.00 km’, RAKE 


JR 29 7] 455.27 km’, , 相 比 2000 年 ( 约 38.19 km?) 


RR tt 


与 塔里木 河 下 游 生态 输 水 量 密切 相关 ;2007 一 2009 年 输 水 量 下 泄 水 量 较 少 ,直接 导致 研究 区 域 在 


2009 年 地 表 水 体 和 季节 性 水 体 的 玫 


积 均 减 少 到 最 低 点 。 
生态 输 水 对 地 表 水 体 变化 的 影响 ,为 区 域 水 生态 保护 提供 科学 依据 。 


研究 结果 有 助 于 全 面 理解 塔里木 河 下 游 
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文章 编号 : 


近年 来 , 受 全 球 气候 变化 和 人 类 活动 的 双重 影 
响 , 塔 里 水 河 从 20 世 纪 50 年 代 开 始 出 现 了 河道 断 
流 、 河 湖水 系 上 肢解 ,湖泊 萎缩 及 水 环境 恶化 等 现象 ， 
危及 到 流域 原本 脆弱 的 生态 环境 ,引发 了 一 系列 生 
态 环境 问题 ,严重 制约 了 区 域 经 济 社会 的 可 持续 发 
feo, ub. A Bere 2001 年 批复 实施 了 《塔里木 
河 近 期 综合 治理 规划 》, 正式 启 动 了 塔里木 河流 域 
综合 治理 工程 。 截 至 2019 年 年 底 , 塔 里 木 河 下 游 输 
水 工程 已 实施 了 20 a, 区域 地 表 水 体面 积 得 到 恢 
复 .地 下 水 埋 深 得 到 明显 抬升 及 生态 环境 出 现 好 转 
等 5 。 其 中 ,以 河 湖 为 主体 的 地 表 水 体面 积 变化 受 
到 塔里木 河 下 游 输 水 工程 的 影响 ,在 参与 塔里木 河 
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下 游 水 循环 .调节 水 平衡 中 扮演 着 重要 角色 ,也 是 
区 域 人 类 最 易 获 取 的 水 资源 。 纵 观 针对 塔里木 河 
下 游 输 水 工程 影响 的 研究 进展 ”“" ,目前 关于 近 20a 
来 塔里木 河 下 游 地 表 水 体面 积 变化 及 其 对 生态 输 
水 响应 的 系统 研究 还 比较 少 。 

2010 年 谷歌 推出 的 Google Earth Engine (GEE) 
云 计算 平台 ,整合 了 多 源 数据 并 结合 了 云 平台 强大 
的 运算 能 力 ,能够 快速 获取 和 处 理 海量 遥感 影像 。 
目前 GEE 已 被 广泛 应 用 于 大 规模 地 理 空间 数据 分 
析 , 为 全 球 和 大 尺度 流域 或 区 域 的 资源 环境 动态 监 
测 、 时 空 变化 特征 解析 提供 了 强大 的 技术 支撑 , 解 
决 了 传统 海量 遥感 影像 处 理 分 析 的 计算 复杂 AS 
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高 昂 等 问题 ”…” ,为 我 们 研究 塔里木 河 下 游 区 域 地 
表 水 体 时 空 变化 提供 了 技术 支持 。 因 此 ,本 文选 取 
塔里木 河 下 游 区 域 为 研究 对 象 , 聚 焦 以 河 湖 为 主体 
的 地 表 水 体面 积 ,采用 GEE 平 台 综合 多 种 统计 分 
析 方 法 ,系统 解析 了 近 20 a 塔里木 河 下 游 区 域 地 表 
水 体面 积 的 变化 特征 ,并 探讨 生态 输 水 对 地 表 水 体 
面积 的 影响 ,为 生态 输 水 后 区 域 水 资源 的 合理 利 
用 .生态 环境 保护 等 提供 科学 依据 。 


1 研究 区 数据 与 方法 


1.1 研究 区 概况 

塔里木 河 下 游 生态 输 水 工程 最 初 是 从 大 西海 
子 水 库 向 下 游 河 道 输 水 , 输 水 量 为 大 西海 子 水 库 的 
下 泄 水 量 ,因此 本 文 的 研究 区 域 为 大 西海 子 水 库 末 
端 (不 包括 大 西海 子 水 库 ) 到 台 特 玛 湖 的 河流 周边 
10 km 的 缓冲 区 域 ( 图 1)。 同 时 ,为 了 准确 解析 生态 
输 水 对 地 表 水 体面 积 时 空 变 化 的 影响 ,将 研究 区 分 
为 上 段 (大 西海 子 水 库 一 英 苏 ,不 包括 大 西海 子 水 


E) .中 段 ( 英 苏 一 阿拉 干 ) 和 下 段 (阿拉 干 一 台 特 玛 
湖 )3 个 区 域 。 
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图 1 塔里木 河 下 游 区 域 概况 


Fig. 1 Overview of the lower reaches of Tarim River 


12 主要 数据 

1.2.1 影像 数据 本 文 主要 基于 GEE 云 计算 平台 ， 
通过 编程 运算 获得 研究 区 域 2000 一 2019 年 所 有 可 
利用 的 Landsat 5, Landsat 7 fll Landsat 821%. 7 m 
本 文 研 究 区 域 的 影像 有 4 景 , 总 计 影像 为 2397 张 ， 


其 中 Landsat 5 为 759 张 、Landsat 7 为 1256 张 和 
Landsat 8 为 582 张 。2012 年 使 用 的 是 Landsat 7 影 
像 数 据 ,2012 年 前 主要 使 用 Landsat 5 和 Landsat 7 影 
像 数 据 ,2012 年 后 主要 使 用 Landsat 7 和 Landsat 8 影 
像 数据 。 

1.2.2 生态 输 水 数据 ”截至 2019 年 年 底 , 向 塔里木 
河 下 游 生 态 输 水 共 20 次 ( 表 1) ,总 计 输 水 81.64x 
10° m ,其 中 塔里木 河 向 下 游 输 水 67.78x10 m C25 i 
输 水 总 量 的 83.03%), 孔 淮河 输 水 13.86 10° m CZ] 
占 16.97%)。 输 水 量 最 大 的 年 份 为 2018 年 ,全 部 由 
塔里木 河 干流 向 下 游 输 水 ,达到 12.15x10’ m ; 输 水 
量 最 小 的 年 份 为 2008 年 , 输 水 量 为 0。 


表 1 2000 一 2019 年 塔里木 河 下 游 生 态 输 水 情况 


Tab.1 Amount of ecological water conveyance in the 


lower reaches of Tarim River from 2000 to 2019 


大 西海 子 水 库 下 泄 水 量 /10' mi 


ii idi am BGA 。 合计 
1 2000 年 5 一 7 月 9923 0 9923 
2 2000 年 11 一 12 月 22655 0 22655 
3 2001 4E4—11 月 34733 3490 38223 
4 20024E.7—11 H 24500 8629 33129 
5 2003 43—11 月 24000 38509 62509 
6 2004 4E 4—6 H 7350 2857 10207 
7 2005 4E 4—11 月 5246 23026 28272 
8 2006 4F.9—11 月 2700 16944 19644 
9 2007 ^F. 10 H 1410 0 1410 
10 2009 4 12 H 1066 0 1066 


11 2010 年 6 一 11 月 0 36393 36393 
12 2011 年 1 一 11 月 0 85211 85211 
13 2012 年 4 一 11 月 0 66716 66716 
14 2013 年 4 一 11 月 0 48769 48769 
15 2014 年 6 月 0 727 727 
16 20154 8—11 H 0 46128 46128 
17 2016 年 8 一 10 月 0 67611 67611 
18 2017 年 4 一 12 月 0 121461 121461 
19 2018 年 2 一 11 月 0 70006 70006 

20 2019 年 8 一 12 月 5002 41351 46353 


13 研究 方法 

1.3.1 地 表 水 体 提取 算法 ”本文 首先 利用 遥感 影像 
数据 计算 研究 区 域 的 归 一 化 植被 指数 (NDVI) JA— 
化 差异 水 体 指 数 (mNDVI) 和 增强 植被 指数 (EVI)， 
计算 公式 如 下 所 示 : 


NDVI= Pnm ~ PRed (1) 
DNIR Prea 


202106.00086v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


Ap 珊 等 : 塔里木 河 下 游 生态 输 水 对 地 表 水 体面 积 变化 的 影响 


mNDWI= P Green PswIR (2) 
P Green + Pswin 
Dun 一 DRed 
EVI= 2.5x (3) 
Pyr + 9-OP Rea — 7-SPBtue + 1 
式 中 > Prine ^ PGreen ^ PRed ^ PNIR 和 PswiR 分 别 代表 蓝 


光波 段 (0.45~0.52 um) 、 绿 光波 段 (0.52~0.60 um) 、 
红 光 波段 (0.63~0.69 jum) 、 近 红外 光波 段 (0.77~ 
0.90 pm) 和 短波 红外 波段 (1.55~1.75 pm)。 

当 mNDWI>EVI 或 mNDWI>NDVI 时 表示 影像 
像素 中 的 植被 信号 比 水 信号 弱 ,说 明 该 像素 可 以 被 
作为 水 体 。 同 时 ,利用 EVI<0.1 这 一 标准 来 尽 可 能 减 
少 湿地 植被 信号 对 湿地 水 体 信 号 提取 的 干扰 "“。 随 
后 ,通过 计算 每 一 幅 遥 感 影像 水 体 出 现 的 频率 来 判 
定 季节 性 水 体 和 永久 水 体 。 对 于 塔里木 河 下 游 的 每 
个 Landsat 像 素 ,用 下 述 公 式 计 算 其 1a 的 水 体 频 率 : 


N, 
FO) = A 3, x 100% (4) 
i=1 


式 中 :下 是 像素 的 水 体 频率 计算 指数 ; y 是 2000 一 
2019 年 的 年 份 ; IN, 是 第 y 年 中 该 像素 的 Landsat 影 
像 所 观测 值 总 数 ; w,, 表示 该 像素 的 一 个 观测 值 是 
否 为 水 体 , 其 中 1 表示 水 体 ,而 0 表示 非 水 体 。 通 过 
计算 水 体 频 率 , 按 照 75% 以 上 为 永久 性 水 体 、25%~ 
75% 为 季节 性 水 体 和 小 于 25% 为 非 水 体 的 规则 六， 
确定 研究 区 域 的 非 水 体 ( 像 元 值 设 为 0) .季节 性 水 
体 ( 值 设 为 1) 和 永久 水 体 ( 值 设 为 2) ,总 的 地 表 水 体 
面积 为 季节 性 水 体面 积 和 永久 水 体面 积 之 和 。 

1.3.2 Mann-Kendall 检验 法 本 文采 用 非 参 数 
Mann-Kendall 趋势 检验 法 (M-K 检 验 法 ) 97 4) gr Hl 
表 水 体面 积 的 时 间 变 化 趋势 。 该 方法 对 待 检验 的 
地 表 水 体面 积 的 时 间 序 列 (TT, T, ,构造 了 一 


个 统计 变量 s rms 


n-l n 
S- >, 5 sign(T, - T) (5) 
i=l jail 
1 ifT,-T,»0 
sign(T, -T)240 if T,- T;-0 (6) 
-1if T,- T,«0 


3X'P : n 29 20(2000—2019 4E) ; sign(x) 为 相应 的 符 
号 函数 。 基 于 “ 待 检 验 地 表 水 体面 积 无 变化 趋势 ” 
的 假设 。 一 般 统计 量 S FE n= 10 时 服从 标准 正 态 
分 布 (均值 为 0) ,方差 由 如 下 公式 计算 
, nn-N2Qn+5)- PDE YQi+5) 
Vila 18 


(7) 


由 此 设计 的 M-K 统 计量 (2 ) 为 


S-l  S>0 
Oo 
Z-40 $20 (8) 
Srl S«0 
Oo 


当 Z>0 时 ,表示 地 表 水 体面 积 呈 增加 趋势 ; 当 
Z<0 时 ,表示 地 表 水 体面 积 呈 减少 趋势 。 当 ZS 
1.96 时 ,表示 通过 w=0.05 的 显著 性 水 平 检验 ; 24 
Z>2.58 时 ,表示 通过 aw= 0.01 的 显著 性 水 平 检验 。 
1.3.3 相关 系数 分 析 ”本 文采 用 皮尔 逊 (Pearson) 相 
关系 数 法 3 于 计算 2000 一 2019 年 塔里木 河 下 游 地 
表 水 体面 积 与 生态 输 水 量 之 间 的 相关 性 。 其 计算 
BRA”: 


n 


245 -Xyv, + Y) 
je (9) 


RP: r 为 水 体面 积 与 输 水 量 之 间 的 相关 系数 ; X, 
和 YY 分 别 为 2 个 变量 的 年 序列 ;全 和 了 为 年 均值 ; n 
为 20(2000 一 2019 年 )。r 处 于 -1 到 1 之 间 , 正 值 表 
示 正 相关 , 负 值 表示 负 相 关 , 其 绝对 值 越 大 ,相关 性 
越 强 , 越 接近 于 0, 相 关 度 越 弱 。 一 般 认为 , [r «0:2 
时 ,为 相关 程度 极 弱 (或 不 相关 ); 0.2 < |r| < 0.4 时 ， 
为 弱 相 关 ; 0.4<|r|<0.6 时 ,为 中 等 程度 相关 ; 0.6 
Ir] < 0.8 时 ,为 强 相关 ; 0.8 < |r| <1 时 ,为 极 强 相关 。 


2 结果 与 分 析 


2.1 地 表 水 体面 积 时 间 变 化 分 析 

图 2 和 表 2 显示 了 2000 一 2019 年 塔里木 河 下 游 
地 表 水 体面 积 变化 的 分 析 结 果 。2000 一 2019 F 
里 木 河 下 游 的 地 表 水 体面 积 .季节 性 水 体面 积 和 永 
久 性 水 体面 积 均 呈现 波动 上 升 趋势 (其 中 地 表 水 体 
面积 变化 具有 显著 性 上 升 趋势 ) ,Z 值 分 别 为 2.37、 
1.65 和 3.08, 上升 的 速率 分 别 为 15.54 km a, 7.12 
km?+ a! $818.41 km?+a!, 20 a 间 地 表 水 体 、 季 节 性 水 
体 和 永久 性 水 体 的 平均 面积 分 别 为 312.21 km’, 
249.05 km 和 63.17 kn? ,最 大 面积 分 别 为 523.93 km? 
(2018 年 )、475.30 km*(2002 年 ) 和 216.44 km?(2017 
年 ) ,最 小 面积 分 别 为 45.89 km?(2009 4E.) 43.38 km? 
(2009 年 ) 和 0.06 km*(2009 年 )。 塔 里 木 河 下 游 区 域 
(大 西海 子 以 下 ) 的 地 表 水 体 .季节 性 水 体 和 永久 性 
水 体 的 面积 从 2000 年 的 49.00 km?, 48.93 km 和 
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— 季节 性 水 体 
(a) 塔里木 河 下 游 区 域 


由 D 
o B 
o o 


地 表 水 体面 积 /km? 
x 
> 
S 


0 1 1 1 1 
2000 2004 2008 2012 2016 2020 
年 份 

507 (c) 中 段 区 域 


上 
= 


地 表 水 体面 积 /km 
N 
© 


10} 

0 ; ; i ， , 

2000 2004 2008 2012 2016 2020 
年 份 


— 永久 水 体 


地 表 水 体面 积 /km? 
上 
e 
© 
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图 2 2000 一 2019 年 塔里木 河 下 游 地 表 水 体面 积 变 化 


Fig.2 Changes of water surface area in the lower reaches of Tarim River from 2000 to 2019 


表 2 2000 一 2019 年 塔里木 河 下 游 地 表 水 体面 积 时 间 变 化 趋势 


Tab.2 Trend of water surface area in the lower reaches of Tarim River from 2000 to 2019 


区 域 类 型 最 大 面积 /km*( 年 份 ) 最 小 面积 km:( 年 份 ) 平均 面积 /km 变化 速率 /km .an Z 值 
整个 区 域 — 季节 性 水 体 475.30(2002) 43.38(2009) 249.05 7.12 1.65 
永久 水 体 216.44(2017) 0.06(2000) 63.17 8.41 3.08" 

地 表 水 体 523.93(2018 ) 45.89(2009) 31221 15.54 2.37 

上 段 区 域 ”季节 性 水 体 17.18(2018) 1.47(2008) 11.21 0.31 243 
永久 水 体 1.08(2017) 0.03(2009 ) 0.23 0.01 0.10 

地 表 水 体 17.65(2018) 1.54(2008) 11.44 0.31 224 

中 段 区 域 。 季节 性 水 体 37.77(2017) 0.20(2008) 19.10 1.58 328" 
永久 水 体 5.44(2018) 0(2000,2002,2009,2010) 0.82 0.16 2.90" 

地 表 水 体 42.50(2018 ) 0.21(2008) 19.92 1.75 3.34” 

下 段 区 域 。 季节 性 水 体 443.15(2002) 29.62(2001) 218.74 5.23 1.59 
永久 水 体 212.14(2017) 0(2000,2001 ) 62.12 8.24 3.05" 

地 表 水 体 463.79(2018) 29.62(2001 ) 280.86 13.48 2.37 


注 :*、** 分 别 表示 在 5%、1% 的 显著 性 水 平 上 显著 。 


0.06 km 增加 到 2019 年 的 498.54 km’ 300.76 km 和 
197.78 km 。 输 水 后 的 第 3 年 (2002 年 ) ,地 表 水 体面 
积 和 季节 性 水 体面 积 出 现 了 急剧 的 增加 ,分 别 从 
2001 年 的 48.58 km? Fil 48.02 km 增加 到 2002 年 的 
491.26 km 和 475.30 km 。2002 一 2006 年 地 表 水 体 
面积 呈现 波动 变化 ,2007 一 2009 年 呈现 急剧 减少 ， 
在 2009 年 地 表 水 体 和 季节 性 水 体 的 面积 均 减 少 到 
最 低 点 ,随后 在 2010 年 出 现 急 剧 增加 ,并 在 2010 一 
2019 年 呈现 缓慢 的 上 升 趋势 。 


塔里木 河 下 游 上 段 的 地 表 水 体 .季节 性 水 体 和 
永久 性 水 体 的 面积 均 呈 现 波动 上 升 趋势 (Z 值 分 别 为 
2.24.2.43 和 0.10) ,上 升 的 速率 分 别 为 0.31 km-a”, 
0.31 km .a 和 0.01 km.a-:。20a 间 地 表 水 体面 积 、 
季节 性 水 体面 积 和 永久 性 水 体面 积 的 平均 面积 
分 别 为 11.44 km 11.21 km? 和 0.23 km’, 最 大 面积 
分 别 为 17.65 km (2018 年 ) .17.18 km (2018 年 fll 
1.08 km? (2017 4E.) ,最 小 面积 分 别 为 1.54 km? (2008 
年 ) 1.47 km*(2008 年 ) 和 0.03 km*(2009 年 )。 与 整 
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个 研究 区 一 样 , 在 2007 一 2009 年 呈现 急剧 减少 ,在 
2008 年 地 表 水 体面 积 和 季节 性 水 体面 积 均 减少 到 
最 低 点 。 

塔里木 河 下 游 中 段 和 下 段 区 域 的 地 表 水 体面 积 、 
季节 性 水 体面 积 和 永久 性 水 体面 积 均 呈 现 波动 的 
上 升 趋势 , 且 除 下 段 区 域 的 季节 性 水 体面 积 外 ,其 
他 水 体 类 型 面积 均 呈 现 显著 性 增加 。20 a 间 中 段 区 
域 地 表 水 体 、 季 节 性 水 体 和 永久 性 水 体 的 面积 增加 
速率 分 别 为 1.75 km…a'、1.58 km2.a 和 0.16 km? a, 
平均 面积 分 别 为 19.92 km?、19.10 km2 和 0.82 kn ,最 
大 面积 分 别 为 42.50 km*(2018 年 )、37.77 km (2017 
年 ) 和 5.44 km' (2018 年 ) ,最 小 面积 分 别 为 0.21 km? 
(2008 年 ) ,0.20 km? (2008 年 ) 和 0 km? (2000, 2002 , 
2009.2010 4E) , 20 a 间 下 段 区 域 地 表 水 体 .季节 
性 水 体 和 永久 性 水 体 的 面积 增加 速率 分 别 为 
13.48 km a1、5.23 km-a” F1 8.24 km*…a ,平均 面积 
分 别 为 280.86 km^,218.74 km 和 62.12 km’, fe H 
积分 别 为 463.79 km?(2018 4E.) 443.15 km'( 2002 4F ) 
fll 212.14 km? (2017 4E) ,最 小 面积 分 别 为 29.62 km? 
(2001 年 )、29.62 km*(2001 年 ) 和 0 km? (2000, 2001 


(a) 2000 年 (b) 2001 年 (c) 2002 年 


1.14% 


(d) 2003 年 


年 )。 与 整个 研究 区 一 样 , 中 段 和 下 段 区 域 的 地 表 
水 体面 积 .季节 性 水 体面 积 在 2007 一 2009 年 呈现 急 
剧 减 少 。 
2.2 水 体面 积 空间 变化 分 析 

塔里木 河 下 游 区 域 地 表 水 体 呈 现 明 显 的 空间 
不 均匀 分 布 , 且 随 着 时 间 的 变化 而 呈现 出 不 同 的 状 
态 (图 3)。 整 体 上 ,水 体 主要 分 布 在 研究 区 域 的 下 
段 。 塔 里 木 河 下 游 区 域 季 节 性 水 体 所 占 的 比重 相 
对 较 大 ,20 a 间 的 平均 占 比 为 83.60% ,最 大 占 比 为 
99.879%(2000 年 ), 最 小 占 比 为 57.31%(2017 年 ); 永 
久 水 体 占 整个 水 体 表面 积 比重 最 大 的 年 份 发 生 在 
2017 4F (42.69% ) , 随后 是 2019 4F (39.67% ) fI 2018 
年 (37.02%)。 图 4 显示 了 塔里木 河 下 游 生态 输 水 起 
始 年 份 2000 年 和 2019 年 的 地 表 水 体面 积 空间 分 
布 ,生态 输 水 引起 的 地 表 水 体面 积 增加 主要 表现 在 
台 特 玛 湖 湖区 (本 文 计算 规 定 台 特 玛 湖 湖区 为 研究 
区 域 中 39°39'54”N 以 下 部 分 )。2000 年 , 台 特 玛 湖 
湖区 只 存在 少量 的 地 表 季 节 性 水 体 ,面积 约 为 
38.19 km2;2019 年 台 特 玛 湖 湖 区 的 地 表 永 久 水 体 和 
季节 性 水 体 均 呈现 显著 的 增加 ,面积 分 别 约 为 


(e) 2004 年 (f) 20054F. 


(g) 2006 年 


(j) 2009 年 


5.45% 


(n) 2013 年 


Si 


5.4396 


8.3196 


(t) 2019 年 图 例 N 
一 河流 人 


永久 水 体 
季节 性 水 体 
其 他 覆盖 类 型 
0 10 km 


注 : 饼 图 代表 季节 性 水 体 与 永久 性 水 体面 积 占 比 。 
图 3 2000 一 2019 年 塔里木 河 下 游 水 体 空间 分 布 


Fig. 3 Annual spatial distributions of water surface area in the lower reaches of Tarim River from 2000 to 2019 


202106.00086v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


tió wA 
267.27 km 和 188.00 km ,总 体 水 面积 约 为 45$.27 km’, 体 在 过 去 20 a 间 的 平均 占 比分 别 为 3.66% .6.38% 和 
相 比 2000 年 增加 了 417.08 kmz( 约 10.92 倍 )。 89.96%( 图 5a) ,说 明 下 段 区 域 拥有 最 大 的 地 表 水 体 
PERN PUR EB .中 段 和 下 段 区 域 的 地 表 水 ”面积 而 上 段 区 域 最 小 ,与 图 3 中 的 空间 分 布 相符 。 
(a) 2000 年 (b) 2019 年 


图 例 
永久 水 体 
季节 性 水 体 
其 他 覆盖 类 型 


图 4 2000 年 和 2019 年 塔里木 河 下 游 水 体 的 空间 分 布 
Fig.4 Spatial distributions of water surface area in the lower reaches of Tarim River in 2000 and 2019 
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(c) 上 段 、 中 段 、 下 段 区 域 地 表 永久 水 体面 积 占 比 (qd) 上段、 中段、 下段 区 域 地 表 季 节 性 水 体面 积 占 比 
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图 5 2000 一 2019 年 塔里木 河 下 游 不 同 水 体 和 不 同 区 域 的 占 比 
Fig. 5 Ratios of water surface area for different types of surface water and different 


regions in the lower reaches of Tarim River from 2000 to 2019 
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其 中 ,20a 间 上 段 .中 段 和 下 段 区 域 的 地 表 永 久 水 体 
的 平均 占 比 分 别 为 0.36% .1.30% 和 98.34% ,季节 性 
水 体 的 平均 占 比分 别 为 4.50% .7.67% 和 87.83% ,说 
明 下 段 区 域 的 永久 性 水 体面 积 在 其 总 的 地 表 水 体 
面积 占 比 中 承担 了 最 大 的 贡献 率 。 网 5$b~d 显 示 了 
2000 一 2019 年 塔里木 河 下 游 上 段 .中 段 和 下 段 区 域 
的 地 表 水 体 .永久 水 体 和 季节 性 水 体 的 比重 变化 。 
整体 上 ,上 段 区 域 的 永久 水 体 和 季节 性 水 体 的 比重 
均 呈 现 减 少 趋势 ,永久 性 水 体 所 占 比重 从 2000 年 的 
100% 下 降 到 2019 年 的 0.36%; 相 反 , 下 段 区 域 的 永 
入水 体 和 季节 性 水 体 的 比重 均 呈 现 增加 趋势 ,永久 
性 水 体 所 占 比 重 从 2000 年 的 0% 上 升 到 2019 年 的 
98.34% , 主要 分 布 在 台 特 玛 湖 湖区 ( 约 占 99.98% ) ; 
中 段 区 域 的 地 表 永 久 水 体 和 季节 性 水 体 的 增长 速 
率 较 低 。 台 特 玛 湖 湖区 水 体面 积 的 扩张 是 下 段 区 
域 水 体面 积 占 比 不 断 增加 的 主要 原因 。 
2.3 地 表 水 体面 积 与 输 水 量 的 相关 性 

近 20a 以 来 ,塔里木 河 下 游 的 地 表 水 体面 积 、 季 
节 性 水 体面 积 和 永久 性 水 体面 积 与 生态 输 水 量 均 


呈现 正 相关 性 ,其 相关 系数 Cr 值 ) 分 别 为 0.45( 正 中 
等 相关 ) .0.15( 正 极 弱 相关 ) 和 0.78( 正 强 相关 ) ,说 
明 塔 里 木 河 下 游 地 表 水 体面 积 ,尤其 是 永久 水 体面 
R, 主要 与 塔里木 河 下 游 输 水 量 相 关 ( 图 6)。 分 区 
域 来 看 ,上 段 .中 段 和 下 段 区 域 的 地 表 水 体面 积 三 
生态 输 水 量 均 呈 现 正 相关 性 ,其 相关 系数 (7 值 ) 分 
别 为 0.68( 正 强 相 关 ) .0.55( 正 中 等 相关 ) 和 0.42( 正 
中 等 相关 ) ;季节 性 水 体面 积 与 生态 输 水 量 也 均 呈 
现 正 相关 性 ,其 相关 系数 (> 值 ) 分 别 为 0.66( 正 强 相 
XH) .0.52( 正 中 等 相关 ) 和 0.09( 正 极 弱 相 关 ) ;永久 
水 体面 积 与 生态 输 水 量 也 均 呈 现 正 相关 性 ,其 相关 
系数 (> 值 ) 分 别 为 0.72( 正 强 相 关 ) .0.62( 正 强 相 关 ) 
和 0.78( 正 强 相 关 )。 这 些 相 关 性 结果 说 明 , 整 体 上 ， 
上 段 区 域 地表 水 体面 积 与 生态 输 水 量 的 相关 性 较 
强 , 而 下 游 区 域 永久 性 水 体面 积 受 生态 输 水 量 的 影 
响 相 对 最 大 。 这 主要 由 于 大 西海 子 下 泄 的 水 量 在 
先后 通过 上 段 和 中 段 区 域 的 河道 后 才能 到 达 下 段 
区 域 的 台 特 玛 湖 湖区 , 且 随 着 输 水 量 的 不 断 增 多 ， 
水 体面 积 不 断 增 大 。 同 时 ,生态 输 水 量 在 2007 一 


(a) 塔里木 河 下 游 地 表 水 体 (b) 上 段 地 表 水 体 (c) 中 段 地 表 水 体 (d) 下 段 地 表 水 体 
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注 : 图 中 红色 直线 是 水 体面 积 与 输 水 量 的 相关 性 拟 合 直 线 , 灰 色 区 域 代 表 拟 合 直 线 95% 的 置信 区 间 。 


图 6 塔里木 河 下 游 区 域 各 类 地 表 水 体面 积 与 生态 输 水 量 的 关系 


Fig.6 Correlations between yearly surface water and amount of ecological water 


conveyance in the lower reaches of Tarim River 
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2009 年 处 于 极 低 值 (其 中 2008 年 输 水 量 为 0) ,直接 
导致 塔里木 河 下 游 区 域 在 2009 年 地 表 水 体面 积 
季节 性 水 体面 积 均 减少 到 最 低 点 。 


3 讨论 


本 文 基于 近 20 a 塔 里 木 河 下 游 的 遥感 影像 资料 
和 生态 输 水 数据 ,解析 了 区 域 地 表 水 体面 积 的 时 空 
变化 特征 ,探讨 了 生态 输 水 对 地 表 水 体面 积 的 影 
响 ,对 深入 理解 塔里木 河 下 游 生 态 输 水 所 带 来 的 生 
态 环境 效应 具有 一 定 参 考 意义 。 塔 里 木 河 下 游 大 
西海 子 水 库 至 台 特 玛 湖 湖区 域 的 降水 稀少 蒸发 
大 ,水 资源 主要 是 大 西海 子 水 库 的 下 泄 水 量 , 即 塔 
里 木 河 下 游 的 生态 输 水 量 。 水 资源 是 干旱 区 生态 
环境 安全 和 经 济 社会 可 持续 发 展 的 首要 限制 性 因 
子 , 随 着 生态 输 水 的 不 断 增加 ,塔里木 河 下 游 地 表 
水 体面 积 的 不 断 扩 大 ,该 区 域 植被 生长 逐渐 得 到 明 
显 恢复 .生态 环境 得 到 持续 改善 “"。 由 于 沿途 河水 
损失 ,使 得 前 几 次 的 生态 输 水 没有 到 达 台 特 玛 湖 ， 
故 台 特 玛 湖 的 永久 水 体 和 季节 性 水 体面 积 在 前 几 
次 输 水 后 没有 显著 增加 ;2007 一 2009 年 大 西海 子 水 
库 的 下 泄 水 量 较 少 ,直接 导致 研究 区 域 的 永久 水 体 
和 季节 性 水 体面 积 均 呈现 显著 下 降 ,并 使 得 地 表 水 
体面 积 在 2009 年 减少 到 最 低 值 ;在 塔里木 河 下 游 上 
段 .中 段 和 下 段 区 域 中 ,以 河 湖 为 主体 的 地 表 水 体 
主要 集中 在 下 段 区 域 的 台 特 玛 湖 湖 区 ,湖区 永久 水 
体 和 季节 性 水 体 的 增加 是 塔里木 河 下 游 水 体面 积 
增 大 的 主要 因素 。 同 时 ,由 于 研究 区 域 的 影像 数据 
可 能 存在 一 定 的 质量 缺陷 和 人 台 特 玛 湖水 体 可 能 受 
到 车 尔 臣 河水 量 的 影响 ,本 文 的 研究 结果 存在 
一 定 的 不 确定 性 ,需要 在 未 来 研究 中 进一步 综合 
虑 其 他 因素 的 影响 。 


4 结论 


通过 以 上 的 结果 分 析 及 讨论 ,本 文 的 主要 结论 
如 下 : 

(1) 2000 一 2019 年 ,塔里木 河 下 游 的 地 表 水 体面 
只 季节 性 水 体面 积 和 永久 性 水 体面 积 均 呈现 波动 
上 升 趋势 ,面积 分 别 从 2000 年 的 49.00 km 48.93 km? 
和 0.06 km 增加 到 2019 年 的 498.54 km^,300.76 km’ 
和 197.78 km2; 输 水 后 的 第 3 年 (2002 年 ) ,地 表 水 体 
面积 和 季节 性 水 体面 积 出 现 了 急剧 增加 ,分别 从 


2001 年 的 48.58 km 和 48.02 km 增加 到 2002 年 的 
491.26 km 和 475.30 kn? ; 上 段 ( 大 西海 子 水 库 一 英 苏 ) 
区 域 的 地 表 水 体 .季节 水 体 和 永久 性 水 体 的 面积 均 
呈现 波动 上 升 趋势 ,增加 的 速率 分 别 为 0.31 km? ar! 
0.31 km-a" 和 0.01 km2.a ;中 有 段 ( 英 苏 一 阿拉 干 ) 区 
域 和 下 段 (阿拉 干 一 台 特 玛 湖 ) 区 域 的 地 表 水 体面 

只 .季节 性 水 体面 积 和 永久 性 水 体面 积 均 呈现 波动 
上 升 趋势 ,其 中 ,中 段 区 域 地 表 水 体 .季节 性 水 体 和 
永久 性 水 体 的 面积 增加 速率 分 别 为 1.75 km? a! 
1.58 km?+ a 410.16 knm.a ,下 段 区 域 地 表 水 体 、 季 
节 性 水 体 和 永久 性 水 体 的 面积 增加 速率 分 别 为 
13.48 km a! .5.23 km'-a  f818.24 km'-a''; 

(2) 塔里木 河 下 游 区 域 的 地 表 水 体 呈 现 出 明显 
的 空间 不 均匀 分 布 ,主要 分 布下 段 区 域 ; 生 态 输 水 
引起 的 地 表 水 体面 积 增 加 主要 表现 在 台 特 玛 湖 湖区 
(39°39'54"N 以 下 部 分 ) ,2019 年 台 特 玛 湖 湖区 的 地 
表 永 久 水 体 和 季节 性 水 体面 积分 别 约 为 267.27 km 
和 188.00 km’, 总 体 水 面积 约 为 455.27 km’, fH EG 
2000 4E ( £4] 38.19 km) 增 加 了 417.08 km? (2% 10.92 
倍 ); 近 20 a OK EB .中 段 和 下 段 区 域 的 地 表 永 久 水 
体面 积 的 平均 占 比 分 别 为 0.36% 1.30% 和 98.34% , 
季节 性 水 体面 积 的 平均 占 比分 别 为 4.50% 、7.67% 和 
87.83%。 

(3) 近 20 a 以 来 塔里木 河 下 游 区 域 的 地 表 水 体 
面积 .季节 性 水 体面 积 和 永久 性 水 体面 积 与 生态 输 
水 量 的 相关 系数 分 别 为 0.45( 正 中 等 相关 ) .0.15( 正 
极 弱 相关 ) 和 0.78( 正 强 相关 ) ,说 明 研 究 区 域 地 表 水 
体面 积 , 尤 其 是 永久 水 体面 积 , 主 要 与 塔里木 河 下 
游 输 水 量 相关 ;其 中 ,中 段 区 域 的 地 表 水 体面 积 与 
生态 输 水 量 的 相关 性 较 强 ,而 下 段 区 域 永久 性 水 体 
面积 受 生 态 输 水 量 的 影响 相对 最 大 ;2007 一 2009 年 
极 低 的 输 水 量 直接 导致 研究 区 域 在 2009 年 地 表 水 
体 和 季节 性 水 体 的 面积 均 减 少 到 最 低 点 。 
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Abstract: Water is a key restricting factor of socioeconomic development and eco-environmental protection in 
the arid inland river basins of northwest China. Because of ecological problems caused by climate change and the 
overexploitation of water resources in the Tarim River Basin (e.g., river breakup, drying-up of lakes, and degrada- 
tion of vegetation cover), a national ecological water conveyance project has been officially implemented in the 
lower reaches of Tarim River since 2000. So far, the region’ s ecological environment has improved steadily, and 
studies have evaluated the impacts (e.g., changes in groundwater level and vegetation cover) of the ecological wa- 
ter conveyance project. However, there are relatively few systematic studies on the changes in the area of surface 
water in the lower reaches of Tarim River in the past 20 years. On the basis of the Landsat 5, 7, and 8 images and 
the annual water conveyance quantity, this paper used the Google Earth Engine, the surface water extraction meth- 
od, and the multiple statistic methods (e.g., the Mann-Kendall trend test and the Pearson correlation coefficient) 
to quantitatively analyze the spatiotemporal changes of the surface water and evaluated the impacts of the ecologi- 
cal water conveyance project on the surface water area. The main results are as follows. Firstly, from 2000 to 
2019, a fluctuating upward trend was detected in the overall surface water area, seasonal water area, and perma- 
nent water area, with increase rates at 15.54 km^-a ', 7.12 km .a ', and 8.41 km .a ', respectively. A fluctuating 
upward trend has also been detected for all three subregions: the upstream region from Daxihaizi Reservoir to 
Yingsu station, the upstream region from Yingsu station to Alagan station, and the downstream region from Alag- 
an station to Lake Taitema. Of these, the downstream region had the largest increase rate for both the seasonal wa- 
ter and permanent water areas, with values at 5.23 km'- a ' and 8.24 km .a |, respectively. In addition, the increase 
in surface water area in Lake Taitema was the main factor for the expansion of surface water area in the lower 
reaches of Tarim River. In 2019, the permanent surface water area, the seasonal water area, and the total surface 
water area of Lake Taitema were about 267.27 km’, 188.00 km’, and 455.27 km’, respectively; the total surface 
water area of Lake Taitema in 2019 was evaluated to have been increased by 417.08 km (about 10.92 times) over 
that in 2000 (about 38.19 km’). Finally, in the past 20 years, there is a strong positive correlation between the sur- 
face water area (especially the permanent water area) and the water quantity of the ecological water conveyance 
project. The water quantity of the ecological water conveyance project was extremely low during 2007—2009, 
which directly caused the lowest values of the total surface water area and the seasonal surface water area. Re- 
sults obtained from this study can contribute to the comprehensive understanding of the impact of ecological wa- 
ter conveyance in the lower reaches of Tarim River and provide a scientific basis for developing regional water 
ecological protection. 

Key words: Tarim River; surface water body; ecological water conveyance; temporal and spatial changes; im- 
pact analysis; Google Earth Engine 


